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Abstract 

In this paper, a family of multi-value cellular automaton (CA) associated with 
traffic flow is presented. The family is obtained by extending the rule-184 CA, which 
is an ultradiscrete analogue to the Burgers equation. CA models in the family show 
both metastable states and stop-and-go waves, which are often observed in real 
traffic flow. Metastable states in the models exist not only on a prominent part of 
a free phase but also in a congested phase. 
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1 Introduction 



Traffic phenomena attract much attention of physicists in recent years. It shows a complex 
phase transition from free to congested state, and many theoretical models have been 
proposed so far[]I]]— 0. Among them we will focus on deterministic cellular automaton 
(CA) models. CA models are simple, flexible, and suitable for computer simulations of 
discrete phenomena. 

The rule-184 CA^] has been widely used as a prototype of deterministic model of 
traffic flow. In the model, lane is single and cars can move by one site at most every 
time step. Its several variations have been proposed recently: First, Fukui and Ishibashi 
(FI) proposed a high speed extension of the rule-184 CA. In the model, cars can move 
more than one site per unit time if there are successive vacant spaces ahead. Second, 
Fuks and Boccara proposed a 'monitored traffic model' 0, which is a kind of quick start 
(QS) model. In the model, drivers can prospect vacant spaces due to car motion in the 
next site and can move more quickly compared with the rule-184 CA. Third, Takayasu 
and Takayasu proposed a slow start (SIS) model |TIJ, in which cars can not move just 



after they stop and wait for a unit time. This represents an asymmetry of stopping and 
starting behavior. We call all these variants a family of rule-184 CA in this paper. 

A desirable condition for CA as a traffic model is that it can show a phase transition 
observed in real data. Let us see an example of real data of Tomei expressway taken by the 
Japan Highway Public Corporation JT2J. Fig{l| shows flow-density diagrams often called 



'fundamental diagram', of both lanes in January 1996 on up line at 170.64 km point from 
Tokyo. Fig.[I| (a) and (b) are diagrams of driving and acceleration lane respectively. We 
see that a phase transition from free to congested state around a density of 25 vehicles/km 
in both lanes, and there is a clearer discontinuity in the acceleration lane. The multiple 
states of flow at the same density around the critical density are also observed particularly 
in the acceleration lane, and a skeleton curve of plot points shows a shape of 'inverse A'. 
These phenomena are often observed in real data|12| . 



There is another observed fact that over a certain critical density, a perturbation to a 
uniform traffic causes a formation of jam, and 'stop-and-go' wave propagates backward fL3| . 



Cars accelerate and decelerate alternately in the wave. We adopt these two experimental 
facts, multiple states and stop-and-go wave, as criteria to judge whether a CA is suitable 
for traffic model or not. Among the rule-184 CA family, only SIS model shows multiple 
states and also shows a stop-and-go wave. 

Recently, a multi- value generalization of the rule-184 CA has been proposed by using 



an ultradiscrete method[14|. Its evolution equation is 



Uj +1 = Uj + min([/' L - Uj) - min(E/f , L - Uj +l ), (1) 



where Uj represents the number of cars at site j and time t, and L is an integer constant. 
Each site is assumed to hold L cars at most. We can prove that if L > and < Uj < L 
for any j, then < Uj +1 < L holds for any j. Thus ([!]) is considered to be (L + 1) -value 
CA. Since (|I|) is obtained from an ultradiscretization of the Burgers equation, we call it 
Burgers CA (BCA)[15|]. Note that BCA is equivalent to the rule-184 CA in a special case 



of L = 1. The positive integer L can be interpreted physically in the following three ways: 
First, the road is an L-lane freeway in a coarse sense, and effect of lane-changing rule is 
not expressed explicitly. In Appendix, we show BCA with L = 2 can be interpreted as 
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a two lane model of the rule-184 CA with an explicit lane-changing rule. Second, Uj/L 
represents a probability distribution of cars in a single-lane freeway. In this case, Uj itself 
no longer represents a real number of cars at site j. Third, also in a single-lane freeway, 
length of each site is assumed to be long enough to contain L cars. By introducing the free 
parameter L into the rule-184 CA family, we can generalize them to multi-value CA and 
obtain rich algebraic properties and wide applications to various transport phenomena. 

Though BCA does not show multiple states at a density in the fundamental diagram, 
we have shown that its high speed extension shows multiple states around a critical 
density (T^j. Metastable states of flow in the model exist on a prominent part of a line 
in a free phase over a critical density. The model is also considered to be a multi-value 
generalization of FI model. This model is mentioned in Sec.|2]3] to compare other models. 

Other models in the rule-184 CA family, SIS model and QS model, are two-value CA 
like FI model. We generalize them to multi-value ones and investigate their properties as 
a traffic model in detail. 

We use a max-plus representation to express time evolution rule of those generalized 
models. This representation comes from the ultradiscrete method which has a close re- 
lation with the max-plus algebra[IT|, as shown in our previous papers. Using max-plus 
representation, we can express evolution rule of the above CA's in a conservation from 
like ([!]), which automatically derives a conservation of total number of cars. 



2 Generalization to multi-value CA 

In this section, we will present four CA models which are multi- value generalization of the 
rule-184 CA family. Before explaining models, we shortly review a rule of car movement 
of BCA. Note that we assume L is a lane number of the road in the description of all 
models. In BCA, cars at site j move to vacant spaces at their next site j + 1 as many as 
possible. Therefore, car flow from j to j + 1 at time t becomes mm(Uj, L — Uj +l ). Thus 
adding min(L r j_ 1 , L — Uj) (flow from j — 1 to j) and — min(L r j, L — Uj +l ) (flow from j to 
j + 1) to Uj (present car number), we obtain C/| +1 and derive (|I|). 

2.1 Multi- value QS model 

First, we generalize QS model to a multi-value one. In the multi-value model, cars at 
site j move to vacant spaces at site j + 1 per unit time as many as possible. This 
movement rule is similar to that of BCA. However, a difference is that drivers at j es- 
timate vacant spaces at j + 1 by predicting how many cars move from j + 1 to j + 2. 
Thus, they estimate vacant spaces at j + 1 at time t to be L — Uj +1 (present spaces) + 
min(L r * +1 , L — Uj +2 ) (predicted spaces due to BCA movement). Therefore considering the 
number of cars coming into and escaping from site j, evolution equation on Uj is given 
by 

U) + min(C/*_ 1 , L-U) + min([/j, L - Uj +1 )) 
- min(E/J, L - U* +1 + min(E/J +1 , L - U] +2 ) ) 

XJ] + min(t/*_ 1 , 2L - U* - L7* +1 ) - min(f/j, 2L - U l J+l - U* +2 ). (2) 



U' +1 = 
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The final form of (0) means that the number of movable cars (t/j) is limited by vacant 
spaces (2L — C/| +1 — Uj +2 ) in their next two sites. In the case of L = 1, this model is 
identical to QS model proposed by Fuks and Boccara, and it is the rule-3212885888 CA 
of a neighborhood size '5' after Wolfram's terminology 0. 



2.2 Multi-value SIS model 

In the original SIS model, lane is single and cars can not move just after they stop and 
wait for a unit time. Movable cars move to vacant spaces in their next site like BCA. 
In the multi-value case, we should distinguish standing cars and moving cars in each 
site because each site can hold plural cars. The number of cars at site j blocked by 
cars at j + 1 at time t — 1 is represented by Uj' 1 — min(L r *~ 1 ,L — Uj+i)- These cars 
cannot move at t, then the maximum number of cars movable to site j at t is given by 
U j ~ {Ut 1 ~ m H U j~\ L - Uj+i)}- Therefore, the multi-value SIS CA is given by 

U* +1 = U t j + min [[/*_! -{U]z\-^m{U]z\,L -U]~ X )},L- uj 

- min [U* - {U]- 1 - min([/j-\ L - U^)}, L - U* +1 )] (3) 

We note that this model includes the original SIS model in the case L — 1. Since U t+1 is 
determined by U l and U t ~ 1 , the evolution equation is second order in time and an effect 
of inertia of cars is included in this model. 



2.3 EBCA2 model 

In our previous paper |l6| , we propose an extended BCA (EBCA) model in which cars 



can move two site forward at a unit time if the successive two sites are not fully occupied. 
In this paper, we call the model 'EBCA2' and another model with a similar extension 
'EBCAl' described in the next subsection. A main difference between two models is a 
priority of movement of fast and slow cars to vacant spaces. Fast cars with speed 2 move 
prior to slow ones with speed 1 in EBCA2 and otherwise in EBCAl. 
Evolution equation of EBCA2 is given by 

U j +1 = U J + min ( & 5-i + 4-2, L - U] + a*„ 1 ) - min(&5 + a)_ x , L - U] +1 + a*), (4) 

where a* = mm(Uj, L — Uj +1 , L — Uj +2 ) which is the number of cars moving by two sites 
and bj = mm(Uj, L — Uj +1 ). In reference [16], we only study the cases of L = 1 and 2. In 
this paper, we show some properties of flow for L in the next section. Note that ([|) in 
the case of L — 1 is FI model which is equivalent to the rule-3436170432 CA. 



2.4 EBCAl model 

There is another possibility when we extend BCA to a high-speed model and we call this 
extended new model EBCAl. In contrast to EBCA2, slow cars with speed 1 move prior 
to fast ones with speed 2. Then car movement from t to t + 1 consists of the following 
two successive procedures, 

a) Cars move to vacant spaces in their next site as many as possible. 
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b) Only cars moved in procedure a) can move more one site and they move to their 
next vacant spaces as many as possible. 

The number of moving cars at site j and time t in procedure a) is given by 6* = min({7 - , L— 
Uj +l ). In procedure b), the number of moving cars at site j + 1 becomes min(6*, L — £/j +2 — 
bj +1 + &j +2 ) ; where the second term in min represents vacant spaces at site j ' + 2 after 
the first procedure a). Therefore, considering a total number of cars entering into and 
escaping from site j, evolution equation of EBCAl is given by 

= r; •//..// 

+ min(6*_ 2 , L — Uj — b)^ + 6$) - min(6*_ 1 , L - U* j+1 - b) + b) +1 ) 
= U^ + mmib^ + b^L-U^ + b^-mm^ + b^L-U^ + b]^). (5) 

Note that (|5|) with L = 1 differs from FI model and it is equivalent to rule-3372206272 
CA. 

We can consider that this rule is an extension of SIS model to high-speed motion since 
cars which stop at procedure a) cannot move at procedure b). Let us consider a relation 
between SIS and EBCAl models in detail. If we write each procedure of EBCAl explicitly 
using an intermediate time step, we obtain 

U J +1/2 = < j ■ b [, i - h 'j ( 6 ) 
U? 1 = U l ; l/2 + min(f^ 1 1/2 - {U^-bU}, L-U? 112 ) 

-min(^ +1/2 - {U* - b)}, L - U%{ /2 ). (7) 

Equations (§) and (^) represent above procedure a) and b) respectively, and Uj +1 ^ 2 denotes 
car number at site j just after procedure a). If we replace t+ 1 by t+ 1/2 in (|l|), we obtain 
Moreover, if we replace t by t + 1/2 and t — 1 by t in (|^), we obtain (g). Therefore 
we can consider EBCAl to be a 'combination' of BCA and SIS rules. 



3 Fundamental diagram and multiple states 

We discuss fundamental diagrams of new CA models described in the previous section. 
In the followings, we will consider a periodic road or a circuit. All models in Sec.|| can be 
expressed in a conservation form such as 

A^ + A jg < = 0, (8) 

where A t and Aj are forward difference operator with respect to indicated variable, and 
represents a traffic flow. Average density p and average flow Q l over all sites are defined 

by 

where K is number of sites in a period. Since all models are in a conservation form, 
average density does not depend on time and we can use p without a script t. 

Figures |2| (a)-(d) are density-flow diagrams of previous models with L = 2 and K = 
30. Figures (a)-(d) corresponds to multi- value QS, multi- value SIS, EBCA2 and EBCAl 
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models respectively. If we set a number N of cars, we obtain a unique density p = N/ KL 
but there are many initial distributions of cars. Each initial distribution does not always 
reach a steady flow and often make a periodic state. Q l can also changes periodically in 
time. Therefore, we plot every points in figures by averaging Q f from t = 2K to t = 4K. 
(We use Q as an average flow.) Note that t = 2K is long enough from initial time to 
obtain a periodic state and 2K time steps is much longer than its period. 

Moreover, we obtain different values of Q from different initial distributions for a given 
p in Figs.|] (b) ~ (d). In Ref.JTB|, this type of state giving different Q values for the same 



density are called 'multiple states' |16[. We use this terminology for other models in this 
paper. Multiple states exist around a critical density. Especially for multi-value SIS model 
(Fig.0 (b)) and EBCAl model (Fig|J (d)), we see a thick branch other than lines forming 
an inverse A shape and the distribution of data in the branch look random. This fact is 
interesting because evolution rules are completely deterministic. 

Next, we focus on EBCAl model and examine its properties in detail. Figure |3| is a 
skeleton diagram of EBCAl model. Branches in Fig.|3| are drawn as straight lines and we 
obtain them using specific initial conditions. For example, states '• • • 11102021110202 • • ■', 
'• ■ ■ 121212 002002 • • ■' and '• ■ ■ 222222 • ■ •' are all steady states in the model, and 

flows of these states are plotted on B, C, D and E respectively in FigfJ. There are two 
branches in congested phase of higher density, that is, B-C and D-E. We call the branch 
D-E and B-C congested phase zero and phase one respectively. It is noted that there are 
many states not on the skeleton in Fig|| For example, state • • -211211 • • • is a steady 
state of the model, and it has density 2/3 and flow 1/2. 

Fig.fl shows a time evolution of flow Q l in EBCAl model with a density of phase 
transition region. Number of total sites K is 60 and 240 in Fig.[| (a) and (b) respectively. 
In the figures, we can see periodic oscillation of flow appearing soon after an initial time. 
We observe that the maximum period of oscillation becomes longer as system size becomes 
larger. Moreover, oscillation in a period is not simple as shown in the figures. Fig.|] (c) 
shows a power spectrum / versus frequency / defined by / = | J2t=o exp(2%ift/T)\/T. 
We set K = 240 and T = 1000. This indicates that the irregular oscillation is similar to 
white noise. 

Fig.^ is a fundamental diagram of EBCAl with L — 1. It is interesting that there exist 
multiple states even in the case L — 1, which is equivalent to the rule-3372206272 CA. 
Among deterministic two-value CA models, only SIS model is known to show metastable 
state so far. Thus, EBCAl with L = 1 becomes the second example of such kind of 
two-value CA. 

At the end of this section, we give a comment on a case of L > 3. Fundamental diagram 
of multi-value QS model is the same regardless of L. In other models, small branches 
increases around the critical density as L becomes larger. Let us show this phenomenon 
using EBCA2 which has clear branches. The fundamental diagram of EBCA2 with L = 7 
is given in Fig|]. There are many branches around the critical density. We can give a 
partial explanation of these branches. Let us assume initial values on all sites are restricted 
to n or L — n where < n < L/2. Then, if we use a transformation V — (U — n)/(L — 2n) 
from U to V, V = corresponds to U = n and V = 1 to U = L — n. We can easily show 
that evolution equation on V is obtained by replacing U by V and L by 1 in (Q). Equation 
(|) with L = 1 is FI model itself and it is two-value CA. Therefore, if all initial values 
of U is restricted to n or L — n, then U at any site is always n or L — n and evolutional 
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state is the same as that of FI model by replacing by n and 1 by L — n. 

Above fact implies a shape of the fundamental diagram. Consider an arbitrary state 
of FI model, that is, (f|) with L = 1. Let us assume density is p and average flow is Q for 
that state. We can obtain a corresponding state of (||) with L > by replacing by n and 
1 by L — n. Then, using (^), we can easily show that density of that state is (1 — ^)p+ t 
and flow is (1 — ^r)Q + ^p. Moreover, we can also show that a fundamental diagram of FI 
model is exactly given by two segments of which end points are (0, 0), (^, |) and (|, |), 
(1, 0) respectively. Therefore, a diagram of EBCA2 with L > has at least a segment 
with end points (f , f ) and (§ + ^, § + ) and that with (| + ^f + |) and (1 - f , f ) 
for any n (0 < n < L/2). Figure ^ is a L = 7 case and the shape of diagram is strictly a 
superposition of all these segments. 



4 Stability of flow 

In this section, we investigate stability of the branches in fundamental diagram obtained 
in the previous section. First, we study a uniform state • - - 11111 ■ • -, which gives the 
maximum flow in each model. Let us define a weak perturbation by a perturbation 
changing a local state 11 to 20, and a strong perturbation by a perturbation changing 1111 
to 2200. Both perturbations clearly do not change the density. These perturbations mean 
that a car or two cars in the uniform flow suddenly decrease their speed and consequently 
fully-occupied sites '2' appear. 

Figures [F] (a)-(c) show an instability of uniform flow by a weak perturbation in multi- 
value SIS, EBCA2 and EBCAl models respectively. We can observe a stop-and-go wave 
propagating backward from the perturbed site in Fig.|7| (a) (multi-value SIS) and (c) 
(EBCAl), while the wave does not appear in (b) (EBCA2). In Fig.[7| (c), flow Q decreases 
by the weak perturbation and transits from A to F in Fig.^. Moreover, we see that the 
final steady state contains locally two states corresponding to B and C in Fig.[| Flow of 
state B is 5/6 and that of C is 1/2, both states are in the phase one, and they give the 
maximum and minimum flow in that phase. Numerical results from various initial states 
with a weak perturbation shows that a state of higher density tends to give this type of 
extreme local states on a branch after some time steps. Moreover, if we give a strong 
perturbation of uniform state • • - 111 • • - corresponding to A, it goes down directly to 
G(l/2,l/2) in the phase zero. Fig.|7| (d) shows an instability due to a strong perturbation 
in EBCAl model. In this case, stop-and-go wave does not occur. Moreover, we see that 
the final state consists of two local states • ■ ■ 200200 • • ■ and ■ ■ ■ 222222 • • ■ corresponding 
to D and E in Fig.|3|, which also give the maximum and minimum flow in the phase zero 
respectively. 

Since a perturbed uniform state becomes a congested one which consists of states 
giving the maximum and the minimum flow in a phase, we can predict a length of the 
final congested bunch. Assume a denotes a ratio of length of final congested bunch to 
length of total sites. Density of B, F and C are 5/12, 1/2 and 3/4 respectively. Thus we 
obtain 

5 M \ 3 1 

— (1 — a) H — a = -, 
12 v ' 4 2' 

and a becomes 1/4 which coincides with that of Fig.|7] (c). In the case of the phase zero, 
we also obtain a = l/4(Fig.[F](d)). 



7 



Contrary to above facts, there also exist some states on branches stable against a 
weak perturbation. Let us consider states on branch A-D in Fig.^]. There exists a state 
• • • 011011011 • • • on the branch and we obtain • • • 011020011 • • • by a weak perturba- 
tion. Against the perturbation, perturbed state is stable and it is still on the branch 
A-D. We call this type of state metastable state. Moreover, let us consider a state 
•• •111120111120- • • corresponding to F in Fig.g We can obtain ••• 201120111120 •• • by 
a weak perturbation, and its flow does not change. Therefore, metastable states exist not 
only on branch A-D but also on the phase one branch. In these cases, the effect of the 
site "2" does not spread over the entire system during the time evolution. Thus some 
part of the multiple branches are metastable states, which can be observed in long time 
simulations starting from random initial conditions. 

5 Concluding discussions 

In this paper, some multi-value generalizations of rule-184 family are studied by using 
max-plus representation. We obtain evolution equations in a conserved form giving (L + 
l)-value CA. We have proved that EBCAl model is a generalization of SIS model, and 
it gives multiple states and stop-and-go waves. We consider that a suitable traffic model 
gives those two phenomena at least. If we observe a real traffic data like Fig.p], existence 
of multiple states is clear. In a real traffic, cars are considered to be always perturbed by 
some traffic effects. Therefore, we consider existence of stable branches against a weak 
perturbation is also desirable for a model. Moreover, stop-and-go waves are often observed 
in a real congested traffic. Considering these points, EBCAl is the most suitable model 
among the models of this paper. 

Next, we give a presumption on a fundamental diagram of real data. In Fig.|I], we 
see the multiple state clearly in an acceleration lane. In a driving lane, since the average 
speed of cars is slower than the acceleration lane, flow tends to be stable and the branch 
is not clear. In the acceleration lane, drivers move faster than cars in the driving lane, 
then over-dense free flow will be likely to occur and we can see the multiple state around 
a critical density in the diagram. Moreover, from the stability analysis we find that there 
are metastable states even in the congested state in the EBCAl model. Although we 
cannot clearly see the fact due to fluctuation of data in Fig.l, we expect that there exists 
both a metastable "weak jam" (line B-C) and a "strong jam" (line D-E) in a real highway 
traffic. 

Finally, we point out future problems. We have investigated mainly about models 
with L = 2. About multiple branches due to larger L, only EBCA2 gives clear branches 
and we can easily discuss their stability. Though other models can also make multiple 
branches, transformation of perturbed states is complicated and its analysis is difficult. 
We should solve this difficulty because we consider EBCAl model is more suitable for a 
traffic model. Moreover, two-dimensional extension of these models and its application 
to a real traffic network is one of the important future problems. 
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Appendix 



In this appendix, we show that BCA model with L = 2 can be interpreted as a two-lane 
model of rule-184 CA. Let us call one of the two lanes A-lane and the other B-lane. Both 
lanes are divided into discrete sites as shown in Fig.Al. Assume Aj and Bj denote number 
of cars at site j and time t in A-lane and B-lane respectively. Since all sites can hold only 
one car at most, A* and Bj are always or 1. A car at site j in A-lane moves according 
to the following rule: 

(a) If site j + 1 in A-lane is empty, the car moves to that site. 

(b) If site j + 1 in A-lane is not empty and site j and j + 1 in B-lane are both empty, 
the car can move to site j + 1 in B-lane. 

(c) Otherwise, the car stays at site j in A-lane. 

As for a car in B-lane, symmetrical rule with respect to lane symbol is applied. The above 
rule can be interpreted as follows: Every car move in its own lane prior to the other lane 
(rule (a)). However, if a car can not move in its own lane and only if there are no traffic in 
the other lane, it changes a lane and move forward in the other lane (rule (b)). Since cars 
move independently in each lane according to rule-184 if rule (b) is omitted, the above 
rule can be interpreted as a two-lane model based on rule-184 with lane-changing effect. 
We can express the rule by a couple of evolution equations as follows: 

A] +1 = A) + min(4„ 1 , 1 - A)) - mm{A), 1 - A) +1 ) 

+ min(l - 4_ x , 1 - A% B)_ x , B)) - min(A$, 1 - B), 1 - (10) 
Bf x = B) + min(5j_ 1 , 1 - B)) - min(Sj, 1 - B] +1 ) 

+ min(l - 1 - B% A)) - mm{B% B] +v 1 - A% 1 - A) +l ). (11) 

The last two terms in both equations express lane-changing effect. If those terms are 
omitted, both equations become independent each other and are equivalent to rule-184 
CA (BCA model with L = 1). If V] is defined by 

E/J = A) + B*, (12) 

it denotes a sum of cars of both lanes. Since A*- and Bj is always or 1, we can derive 
an evolution equation on Uj from ( |10|) and (|TT|) as follows: 

U] +l = U) + min([/|_ 1 , 2 - E/J) - min(C/* ) 2 - £/J +1 ). (13) 
This equation is equivalent to BCA with L = 2. 
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Figure Captions 

Fig.l Observed data of flow (vehicles/5min.) versus density (vehicles/Km) on Tomei 
expressway. This diagram was taken by Japan Highway Public Corporation, (a) 
driving lane, (b) acceleration lane. 

Fig. 2 Fundamental diagrams of multi-value CA models with L = 2 and K = 30: (a) 
multi-value QS, (b) multi-value SIS, (c) EBCA2, (d)EBCAl models. 

Fig. 3 Schematic fundamental diagram of EBCAl model. 

Fig. 4 Time evolution of flow in EBCAl model, (a) K = 60, (b) K = 240, (c) power 
spectrum I versus frequency / of traffic flow for K = 240. 

Fig. 5 Fundamental diagram of EBCAl model with L — 1. 

Fig. 6 Fundamental diagram of EBCA2 model with L = 7. 

Fig. 7 Instability of uniform flow due to a weak perturbation, (a) multi- value SIS, (b) 
EBCA2, (c) EBCAl models, (d) shows an instability due to a strong perturbation 
in EBCAl model. In all figures, black, dark gray, light gray squares denote a value 
2, 1, respectively. 

Fig. Al Interpretation of BCA with L = 2 as a model of coupled single lanes. 
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Figure 1 by K. Nishinari and D. Takahashi 
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Figure Al by K. Nishinari and D. Takahashi 
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Figure 2(a) by K. Nishinari and D. Takahashi 
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Figure 2(b) by K. Nishinari and D. Takahashi 




Figure 2(c) by K. Nishinari and D. Takahashi 
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Figure 2(d) by K. Nishinari and D. Takahashi 
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Figure 3 by K. Nishinari and D. Takahashi 




Figure 4(a) by K. Nishinari and D. Takahashi 




Figure 4(b) by K. Nishinari and D. Takahashi 




Figure 4(c) by K. Nishinari and D. Takahashi 




Figure 5 by K. Nishinari and D. Takahashi 
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Figure 6 by K. Nishinari and D. Takahashi 
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Figure 7 by K. Nishinari and D. Takahashi 




